Bright and dark solitary waves in a one-dimensional spin-polarized gas of fermionic 
atoms with p-wave interactions in a hard-wall trap 
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In this paper we elucidate the physics underlying the fact that both bright and dark solitary 
waves can arise in a one-dimensional spin-polarized gas of fermionic atoms with attractive three- 
dimensional p-wave interactions in a hard-wall trap. This is possible since the one-dimensional 
fermion system can be mapped to a system of bosons described by the Lieb-Linger model with 
either repulsive or attractive delta-function interactions which can support solitary waves. 

PACS numbers: 03.75.-b,05.30.Fk 



Advances in experimental techniques for probing ul- 
tracold gases have resulted in a shift in emphasis in the- 
oretical and experimental work in recent years from ef- 
fective field approaches to more refined methods capable 
of dealing with strong correlations. Such strong corre- 
lations occur in ultracold gases confined in de Broglie 
waveguides with transverse trapping so tight that the 
atomic dynamics is essentially one-dimensional (ID) 
with confinement-induced resonances P, Q allowing Fes- 
hbach resonance tuning Q of the effective ID interactions 
to very large values. This has led to experimental verifi- 
cation [4|, i6,] of the fermionization of bosonic ultracold 
vapors in such geometries predicted by the Fermi-Bose 
(FB) mapping method Q, an exact mapping of a ID gas 
of bosons with point hard core repulsions, the "Tonks- 
Girardeau" (TG) gas, to an ideal spin-aligned Fermi gas. 
The "fermionic Tonks-Girardeau" (FTG) gas [l,!], a ID 
spin-aligned Fermi gas with very strong attractive inter- 
actions, can be realized by a 3D p-wave Feshbach reso- 
nance as, e.g^in ultracold '^''K vapor [loj. It has been 
pointed out [1, @, |^ that the generalized FB mapping 
0, B B [m can be exploited to map the ideal FTG gas 
with infinitely strong attractive interactions to the ideal 
Base gas, leading to "bosonization" of many properties 
of this Fermi system. In a recent paper l^] we examined 
how the properties of such an FTG gas on a mesoscopic 
ring are changed for an even number N of fcrmions when 
the strength of the atom-atom attraction is made finite, 
in which case there exists an FB mapping between the 
FTG gas and a system of N bosons with repulsive delta- 
function interactions. 



Liniger (LL) model and spin-aligned fermions with 
attractive interactions of a generalized FTG form with 
reciprocal fermionic coupling constant gf^; here sgn(x) 
is -1-1 (— 1) if x > (a; < 0). All previous work, includ- 
ing 12] and a recent Bethe ansatz treatment of such a 
system in a hard- wall trap [13] , have considered only the 
case where the strongly attractive spin-aligned Fermi gas 
maps to the LL model with repulsive interactions, corre- 
sponding to a negative ID scattering length aio- This 
restriction is not necessary and a minor change in the 
parameters results in a mapping of the FTG gas to an 
LL model with attractive interactions, corresponding to 
a positive scattering length < anj < oo, an equally 
physical regime. The ideal FTG gas |aiD| — s- oo lies 
on the "knife edge" between the attractive and repul- 
sive LL regimes, where am jumps discontinuously from 
— oo to -|-oo, as might have been expected from the fact 
that it corresponds to a zero-energy scattering resonance. 
The purpose of the present paper is to elucidate how this 
works, and to compare and contrast the ground states 
of the spin-aligned Fermi gas in these two regimes. In 
particular, for the case of a hard-wall trap and using the 
mapping between the FTG and LL model, we expose that 
both bright and dark solitary waves can arise. The possi- 
bility of solitons in degenerate Fermi gases has 
certainly been investigated previously, and the novelty in 
the present work is that it includes the p-wave fermion- 
fermion interactions, highlights the relation between the 
Fermi and Bose systems, and exposes the role played by 
the solitary waves of the underlying Bose system. 



It has been realized in recent years that the FB map- 
ping method used to exactly solve the TG gas and 
FTG gas d, 0] is not restricted to the TG case of 
point hard core boson-boson repulsion and the ideal FTG 
case of infinite zero-range fermion-fermion attraction. In 
fact, the same unit antisymmetric mapping function used 
there, A{xi,- ■ ■ ,xn) = ni<j<£<Ar sgn(a;j -Xf), provides 
an exact mapping between bosons with delta- function 
repulsions gijj6{xj — xi) of any strength, i.e., the Lieb- 



3D p-wave resonance and induced ID interaction: 3D 
p-wave Feshbach resonances have been observed in some 
species of ultracold gases, e.g., in [l^. When such 
a gas is contained in a de Broglie waveguide with tight 
transverse trapping so as to reach the regime of effectively 
ID dynamics, such a 3D resonance leads to a nearby 
confinement-induced resonance (CIR) in the ID scatter- 
ing, and it has been shown [1, H, [9| that the 3D p-wave 
scattering length ap and ID scattering length oid in the 
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neighborhood of such a resonance have the connection 



6v; 



a,D = ^[l + 12(Fp/ai)|C(-l/2,l) 



(1) 



^y h/mrcd^± is the transverse oscillator 
ip = — limfc^o tan5p(fc)/fc'^ is the p-wave 



where a± 
length, Vp 

"scattering volume" , ap is the p-wave scattering length, 
C(-l/2,l) = -C(3/2)/47r = -0.2079... is the Hurwitz 
zeta function evaluated at (—1/2, 1), and rrired = m/2 is 
the reduced mass. The expression ([T]) has a resonance 
at a negative critical value V^"* /a\ = —0.4009- • •, im- 
plying that the CIR only occurs when ap < 0. However, 
by varying an external magnetic field so as to sweep ap 
through the critical value, one can in principle vary gid 
from —oo to -t-oo. The ideal FTG gas, which maps to the 
ideal Bose gas, correspond to the point where aiD jumps 
discontinuously from — oo to +oo, and nearby values of 
Op correspond to an "imperfect FTG gas" , which we shall 
call herein simply an FTG gas. am is defined implicitly 
by the contact condition [3, M 0| 



*F(a;i2 = 0+) = -*F(a;i2 = 0-) = -aiD*i.(a;i2 ^ 0±) 

(2) 

where '^f{xi2) is the fermionic relative wave function 
and the prime denotes differentiation. 

Several different ID pseudopotentials have been used 
in the literature to represent a ID potential leading to 
this contact condition, but most of these involve com- 
binations of delta functions and derivatives which make 
them physically opaque. However, the physics is clearly 
exhibited if one instead uses a suitable zero-range limit 
of a deep and narrow square well of depth Vq and width 
2xq, with the zero-range limit Vq ^ oo and — > car- 
ried out at constant VqXq d, [^, Consider first the 
two-body problem A'' = 2 in infinite space. The relative 
wave function ^f{xi2) of the ground state inside the well 
is sin(KXi2) where the value of the positive parameter k is 
determined by the total energy (potential plus kinetic), or 
equivalently, by aiD- One finds kxq — as ^ 

[8|, 19[. The ideal FTG gas (total energy zero, i.e., zero- 
energy resonance) corresponds to infinite oinand hence 
KXo = f . The imperfect FTG gas studied in [1^ has neg- 
ative aiD and positive total energy, i.e., inside the well 
the positive kinetic energy exceeds the negative potential 
energy — Vq by a finite amount. What we wish to point 
out here is that exactly the same relation kxq = f + ^aio 
also applies to the case where aio is positive and the neg- 
ative potential energy exceeds the positive kinetic energy 
inside the well, resulting in a finite and negative total 
energy, i.e^ a bound ground state. The generalized FB 
mapping [3, H, [111 *B(a;i2) = sgn(xi2)^'F(xi2) maps 
the fermionic wave function to a bosonic one with 
the same scattering length oid. The case aio < of this 
mapping was illustrated in Fig. 1 of [I4I, and the case 
diD > which we wish to study here is illustrated by the 
present Fig[T] For consistency of the mapping one should 
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VB(Xi2)=giD5(Xi2) 

iB=sgn(Xi2)ilJp 
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ilJp= sgn(xi2)ilJB 



FIG. 1: Two-particle fermionic relative wave function "ifp 
and mapped bosonic wave function 'I's as a function of 
^12 = X1—X2, for the case where the negative potential energy 
slightly exceeds the positive kinetic energy, so the strongly at- 
tractive Fermi system maps to a Lieb-Liniger Bose gas with 
weak attractions. The potential for both 'I/f and ^l/s is a 
deep and narrow square well plus a point hard core, but in 
the zero-range limit its effect on outside the well is the 
same as that of Vb ~ g^o^{^i2), where here g^jj < 0. 



add a zero-diameter hard core at X12 = to the poten- 
tial, which has no effect on "^p In the zero-range 
limit the effect of this potential on '^b outside the well 
is exactly the same as that of an LL delta function in- 
teraction Vb = gii:)S{xi2) with g^j^ — —2h^/maiD < 0. 
The corresponding fermionic coupling constant is recip- 
rocally related to g^, i.e., g[^ = —2h'^aiB)/m < and 

9id9id ~ 7^ fl^- Outside the well the fermionic and 
mapped bosonic wave functions are both decaying expo- 
nentials, i.e., the ground state is bound. In the case of 
enclosure in a longitudinal box or well, or trapping on a 
ring, the wave functions are modified at large distances 
in accordance with the boundary conditions and/or trap 
potential, but the above represents the exact behavior of 
the wave functions as X12 — > 0. 

Many-fermion ground state: For a system of N 
fcrmions the Fermi wave function may be related to the 
underlying Bose wave function via the mapping "^p = 
A'i'B, where is the ground state of the N-boson 
LL Hamiltonian. The relative strengths of the induced 
ID fermion interactions and the LL model interactions 
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are quantified in the dimensionless coupling coefRcients 
If = "^T" and 75 = , respectively, where n 



N 
L 

The ideal 



is the linear atomic density, and jfJb = 
FTG gas is realized in the limit aio —00 so that 
7;^ — > 00 and 75 0. Here we are interested in the limit 
of finite but very large |ai_D|, implying \^p\ >> 1 and 
I 7b I << 1. This is precisely the limit where the ground 
state properties of the LL Hamiltonian can be accurately 
captured using mean field theory in which all N bosons 
are assumed to occupy the same normalized single par- 
ticle orbital ip that is determined as the ground state 
solution of the Gross-Pitaevskii equation (GPE) [ijj 



2m dx'^ 



-5fz,(iV~l)|0p0, 



(3) 



with /i the chemical potential, and the GPE is to be 
solved subject to the hard-wall boundary conditions 
0(0) = 4>{L) = for a trap of length L. (We employ 
a hard-wall trap due to the fact that the solutions of 
Eq. ([3]) are well known, but our general findings would 
also apply to a harmonic trap). The iV-particle fermion 
wave function can then be written as 4'i?(xi, • • • , xn) = 

A{xi, ■ ■ ■ , XN)Y[j=i 'f'i^j)^ ^^'^ corresponding re- 
duced one-body density matrix (OBDM) may be ex- 
pressed as pi{x,x') = N(f>{x)(t)* {x')[F{x,x')]^~^ where 
F{x,x') = J^^sgnix - x")sgnix' - x")\^{x")\^dx" 
The single-particle momentum spectrum n{k) for the 
Fermi gas may be obtained by Fourier transformation 
of the OBDM and whereas the density p{x) = 

pi{x,x) = A^|0(a;)p is the same for both the Fermi and 
mapped Bose systems, the momentum spectra of the two 
systems differ greatly [14 1. 

fliD < 0; In this case gf^ > and the mapped LL bo- 
son interactions are repulsive, 5^ > The ground 
state solution of the GPE ^ with hard-wall boundary 
conditions is a dark solitary wave whose spatial form is 
an elliptic Jacobi snoidal function: Exhaustive details 
of the form of the solutions are given in Ref. 2l|. In 



the limit of a single atom the GPE ground state de- 
generates into a half period of a sine function over the 
length L of the trap. In contrast, for a large number of 
atoms > > 1 the solution tends to the value y/n within 
the center of the trap, but approaches the classic hyper- 
bolic tangent form of a dark soliton near the hard wall 
boundaries within a distance characterized by the heal- 
ing length = •\/|ai£)|/(2n) << L which follows from 
the relation h^/(2m£,l) = gf^n [21]. 

Figure [2] shows the calculated scaled density profile 
p{x)L for = 14 fermions and 73 =0.1 for the case 
a-iD < (dashed line) and also the ideal FTG gas (sohd 
line) and we see that the density profile is broadened 
by the finite fermion interactions. For the ideal FTG 
gas the competing effects of the infinitely strong attrac- 
tive ID interactions and repulsion due to Fermi degen- 
eracy pressure exactly cancel and that is why the ideal 
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FIG. 2: Scaled density profile p{x)L versus scaled coordinate 
x/L for the ideal FTG (solid line), finite aio < (dashed 
line), and finite am > (dotted line). The results are for 
iV = 14 and |7b| = 0.1. 



FTG maps to the ideal Bose gas with no residual interac- 
tions. In contrast, for the case considered here with finite 
but large am < 0, the Fermi degeneracy pressure domi- 
nates over the attractive interactions and that is why the 
FTG is mapped to a LL model with repulsive interac- 
tions that broaden the density profile. This broadening 
of the density profile in the presence of decreasing attrac- 
tive induced ID fermion interactions was first found by 
Hao, Zhang, and Chen lj| using exact solutions for the 
fermion system based on the Bethe ansatz. Here we have 
elucidated the underlying physics by exposing the rela- 
tion to the mapped LL model with repulsive interactions 
and its dark solitary wave solution. 

Figure [3] shows the scaled momentum spectrum 
27Tn{k)/L versus scaled momentum kL/2'K for \^b\ = 0.1, 
= 14 fermions, and am < (solid line). Bender, 
Erker, and Granger (20| have previously shown that for 
the ideal FTG the high momentum tails of the momen- 
tum spectrum are intimately related to the decay of the 
OBDM pi{x,x') as one moves off-diagonal x ^ x' . In 
the thermodynamic limit they found that the decay of 
the correlations is exponential with decay constant 2n, 
and in our notation the scaled momentum spectrum be- 
comes 



27rn(fc) 



{Am? 



(4) 



Thus, the larger the density n is the more rapidly the 
correlations decay, and the larger the width Afc = 4n of 
the momentum spectrum is. The dotted line in Fig. [3] 
shows the momentum spectrum using Eq. (j4]) and yields 
reasonable agreement with the high momentum tails in 
the numerical solution (solid line) even though we have 
a finite trap and are not in the thermodynamic limit. 
o-iD > 0.' In this case < and the mapped LL bo- 



4 



Irrr 




5 10 15 20 25 

WU2n 



FIG. 3: Scaled momentum spectrum 27rn(fc)/L versus scaled 
momentum kL/2n for |7fl| = 0.1, A'^ = 14 fermions, am < 
(solid line), am > (dashed line), and using Eq. Q (dotted 
line). 

son interactions are attractive, < [l^]. The ground 
state solution of the GPE ([3]) with hard-wall boundary 
conditions is a bright solitary wave whose spatial form is 
an elliptic Jacobi cnoidal function: Details of the form 
of the solutions are given in Ref. jS^]. In the limit of 
a single atom the ground state solution degenerates into 
a half period of a sine function over the length L of the 
trap, whereas for a large number of atoms >> 1 the 
solution becomes a localized peak at the center of the 
trap and approaches the classic hyperbolic secant form 
of a bright soliton with a characteristic width w < L de- 
termined by h^/{2mw^) = \gfjj\N/w Thus, for a 
large number of atoms the bright soliton width becomes 
independent of the trap length L, w = aiD/(4iV). 

Figure [2] shows the calculated scaled density profile 
p{x)L for iV = 14 fermions and \^b\ =0.1 for finite 
0.1D > (dotted line) and also the ideal FTG gas (solid 
line) and we see that the density profile is narrowed. This 
can be intuited using the fact that the mapped LL model 
has attractive interactions so that the density profile is 
narrowed compared to the free Bose gas, which maps to 
the ideal FTG. 

Figure [3] shows the scaled momentum spectrum 
2TTn{k)/L versus momentum kL/2'K for \^b\ = 0.1, N = 
14 fermions, and aiD > (dashed line), and by compar- 
ison with the case auj < (solid line) the momentum 
spectrum is broader. This may be understood as fol- 
lows: Inspection of the plots (not shown) of the OBDM 
pi{x,x') for aiD > shows that it decays significantly 
more rapidly as one moves off-diagonal x x' in com- 
parison to the case aiu < 0, which implies that the mo- 
mentum spectrum will be broader [2(]| , see the discussion 
surrounding Eq. This is also consistent with the fact 
that the peak density of the bright solitary wave is higher 
than that for the dark solitary wave, see Fig. [21 and the 
momentum spectrum width Afc increases with density. 



In summary, we have shown that both bright and 
dark solitary waves can arise in a one-dimensional spin- 
polarized gas of fermionic atoms with attractive 3D p- 
wave interactions in a hard-wall trap, by virtue of the 
fact that the fermion system can be mapped to a sys- 
tem of interacting bosons with induced ID interactions 
that can be either attractive or repulsive, exposing the 
underlying solitary waves. 
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